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ABSTRACT 

To assess the resolution of the local gravity field from kinematic measurements, a state 
model for motion in the gravity field of the earth is formulated. The resulting set of equations can 
accommodate gravity gradients, specific force, acceleration, velocity and position as input data and 
can take into account approximation errors as well as sensor errors. 

1. PROBLEM STATEMENT 

The last few years have seen major advances in kinematic methods of gravimetry. 
Shipbome gravimetry, already a reliable tool, will be further enhanced by using accurate position 
and velocity information from differential GPS. Airborne gravimetry in either the fixed wing 
aircraft or helicopter mode experienced a resurgence over the last few years and is at the point 
where it provides gravity information of acceptable accuracy for wavelengths down to 10 or 15 
km Airborne gravity gradiometry has entered the testing stage and holds great potential for short 
wavelength resolution. Compared to even fifteen years ago, there is now a variety of sensors on 
the market and it appears that a judicious combination will yield information on different parts o 
the gravity spectrum. To assess different sensor configurations, a model is needed which allows 
the combination of kinematic measurements from gravity gradiometers, dynamic gravity meters, 
inertial sensors, differential GPS, laser altimeters, precise pressure altimeters and similar devices. 
The model must allow for the interaction of gravitational and inertial measurements and must be 
able to take sensor biases and measurement noise into account. The formulation of such a model 
using state space techniques is the topic of this extended abstract. A detailed derivation with a 
comprehensive list of references will be published in the near future. 

2. THE STATE SPACE MODEL OF KINEMATIC GEODESY 

Newton's second law for motion in the gravitational field of the Earth, expressed in an 
inertial frame of reference (i), will be taken as the starting point 


= f i + g i 


( 1 ) 


where r; is the position vector from the origin of the inertial frame to the moving object and is 
the second time derivative of this vector, f, is the specific force vector, and g i is the vector of all 

gravitational accelerations acting on the moving object. , . f 

^ The set of nonlinear second-order differential equations (1) can be transformed into a set of 
first-order equations of the form 
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In general, measurements will not be taken in an inertial frame of reference but in an arbitrary body 
frame (b). They can, be transformed into an inertial reference frame by 

fi = Rib fb (3) 

where Rih is a three-dimensional orthogonal matrix transforming f b from the body frame (b) to 
the inertial frame (i). Note that the subscripts denote the direction of the rotation, not the element in 
the matrix. It is obvious from equation (3) that measurements f b can only be used if Ri b is known 
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or can be measured. Thus, a set of three first-order differential equations for the rotation rates R ih 
has to be added to equation (2). They are of the form b 


Rib ~ Rib G 


ib 


where £2^ is the skew-symmetric matrix of angular velocities. 

t___ Similarly, g ?' aV u itati , < ? n L is J usuall y not Sf ven in an inertial frame but in the Conventional 
Terrestrial frame, which will be denoted by (e) in the following. We thus have the transformation 

8 i = R ie g e ♦ ( 5 ) 

n^Uv rot f tlon rate of the Earth can be considered constant for the applications discussed here, 

smte^uations eqUatl0nS ^ R ' C h3Ve t0 be added ' Usin 8 equations (3) to (5) in (2) leads to the 


(4) 


( 6 ) 


^anshdonar parameter.^ 101 ' 011 “ three ' dimensional Euclidean space by three rotational and three 

The implicit assumption in equation (6) is that gravitation g e is known. In that case the 
measurement of the specific force f b and the rotation rates Rj b is sufficient to determine position 
velocity and attitude as functions of time. If gravitation is not known as a function of time, 
ltional measurements are necessary. They can be of two types, inertial or gravitational. In the 
first case, v, is determined independently, in the second case g j is obtained. In both cases 
gravitation and inertia can be separated in the second set of equations in (6). Thus in principle’ 
position, velocity, attitude and gravitation can be obtained by measuring f b , R ib and either v or e 
independently. The second case which is that of gravity gradiometry will now be discussed. ' 
f ° tra . nsfo ^ n e< Juation (6) into a system that also admits gravity gradiometer measurements, 

a , s . et f .? ua !'? f ns Scribing the change of the gravitational vector with time has to be added It is 
obtained by differentiating equation (5) 
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g i - R ie g c + Rie g e = Rie (Qjf g e + g e ) 
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is the matrix of second-order gradients of the gravitational potential V e 
includes gradient measurements of type (7) is therefore of the form 


The state vector which 
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3 CHANGE OF COORDINATE FRAME AND LINEARIZATION 

In many cases results are required in an earth-fixed coordinate frame and transformation of 

coiTespondfog states'in^a r^^ingTram^^The ^^ V6n 

here without proof for the conventional terrestrial frame (e) and for the local level tr t )■ 

In the first case, we get 
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where g e is the gravity vector. The transformation into the local-level frame 1, using curvilinear 
(ip, X, h)-coordinates, follows along similar lines, and results in the state vector 


*1 = 


where 


D = 


n 

VI 
Rib 
^ 81 

1/R 

0 

0 


Dvj 

R lb f b - M® + Of 1 ) V! + g, 

Rib of 


(ID 


0 

1/(R cos cp) 

0 


0 

0 

1 


The state equations derived up to this point are all nonlinear. The first step towards the 
solution of the differential equation system is usually 

reference trajectory is obtained by introducing a gravity model into equation (9) j and by integrati g 
the gravity corrected measurements. Equation (9) can be written in the general form 

x = f {r, v, w, g ) (12) 

where w denotes the angular velocities in the Q* and the & matrices. Subscripts have been 
omitted for convenience. It can be rewritten as 

io + dx° = f (r° + dr, v° + dv, w° + dw, g° + dg } (13) 

where the superscript o denotes the reference trajectory and d the perturbation. By separating the 

reference trajectory 

x° = f {r°, v°, w°, g° ) 

from equation (12) and considering only first-order perturbations, the following set of equations is 
obtained 
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fdr\ f Fu F 12 F 13 F m N /<lr\ 

dV _ F21 - - - dv 

dw " F 3 , - - - dw ( 15 ) 

<dgj V F 4 , -- — F44 J ydg^ 

which is a set of linear homogeneous state equations of the form 

x = Fx . ( 16 ) 

The matrices Fn to F33 are well known from the literature. F41 to F44 are obtained by developing 
the perturbation solution for g. In the local-level frame, we get the reference solution 

g° = {G° (r°) - Q' e (r°) Q>e( r o)} v ° + Qe!( r o) gO(,-o) ( 17 ) 

and the perturbation solution 

dg = {G° - O lc O ie } dv + dG v° - O cl dg - dQ el g° ( 18 ) 


where 


nr tr ro\ j. f dG° -■ . dG° dG° \ 

dG _ (G - G°) + (-jp-dr +-^-dv +-j^-dwj 


Using equations ( 18 ) and ( 19 ) in ( 15 ) will give the desired solution. 


SENSOR ERROR MODEL 


The linearized model given by equation ( 16 ) is the kinematic description of rigid body 
motion in the gravity field of the Earth. Except for approximation errors, it is a rigorous model. If 
sensors are used to determine the trajectory and the gravity vector, the above model has to be 
augmented by a set of error states for each of the sensors used and by system noise. The resulting 
equations will be of the general form 



where xj contains the states in equation ( 15 ) and x 2 the error states, and where Fj and F 2 are the 
corresponding dynamics matrices, and Gu describes the system noise. The F 2 -matrix has a block- 
diagonal structure because the error states of the individual sensors are uncorrelated. 

The error states are modeled in a stochastic manner by random biases, random walk 
models, or simple Gauss-Markov processes. Some judgement is required to decide whether a 
specific error source is modeled into x 2 or into Gu. If real-time estimation is performed with the 
model (20), a small state vector is often desirable. 
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